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EXECUTIVE  SUMMARY 


THE  MAGNETIC  INDUCTION  OF  THE  SYSTEM  CONSISTING  OF 
A  COIL  AND  A  FERROMAGNETIC  PROLATE  SPHEROIDAL  SHELL 

OBJECTIVE 

The  objective  of  this  work  was  to  derive  solutions  to  static  ferro¬ 
magnetic  problems  that  included  both  current-carrying  coils  and  linear 
ferromagnetic  bodies.  The  solutions  are  intended  for  comparison  with 
solutions  to  ferromagnetic  problems  obtained  by  various  americal  tech¬ 
niques  such  as  the  finite  difference  method,  the  finite  element  method, 
and  the  integral  equation  iterative  solution  method. 

APPROACH 

After  deriving  the  governing  differential  equation  from  Maxwell's  . 
equations  for  classical  magnetostatic  field  theory,  the  method  of  separa¬ 
tion  of  variables  was  employed  to  obtain  the  problem  solution. 

RESULTS 

The  magnetic  induction  was  calculated  for  two  geometries  (configura¬ 
tions)  of  a  ferromagnetic  prolate  spheroidal  shell  and  a  current-carrying 
conductor.  The  first  case  was  for  an  infinitesimally  thin  current  band 
carrying  a  stationary  current  and  surrounding  the  spheroidal  shell.  The 
second  case  was  for  a  spheroidal  shell  surrounding  an  infinitesimally  thin 
current  band.  The  ferromagnetic  bodies  were  assumed  to  be  linear  and 
homogeneous.  The  reduction  of  the  solutions  to  that  of  a  current  band  in 
free  space  is  shown  when  the  permeability  of  the  ferromagnetic  spheroidal 
shell  is  allowed  to  approach  that  of  free  space. 

RECOMMENDATIONS 

It  is  recommended  that  the  derived  solutions  be  programmed  on  a  digi¬ 
tal  computer  for  direct  comparison  of  these  results  to  those  obtained  by 
various  numerical  methods.  There  are  plans  to  implement  these  recommenda¬ 
tions  during  the  fiscal  years  1979  and  1980. 
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ABSTRACT 


The  magnetic  induction  is  calculated  for  two 
configurations  consisting  of:  (1)  an  infinitesi¬ 
mally  thin  prolate  spheroidal  current  band 
carrying  a  stationary  current  surrounding  a  fer¬ 
romagnetic  prolate  spheroidal  shell,  and  (2)  an 
infinitesimally  thin  band  internal  to  a  ferro¬ 
magnetic  prolate  spheroidal  shell.  The  ferro¬ 
magnetic  body  is  assumed  to  be  linear  and 
homogeneous.  The  reduction  of  the  solutions  to 
that  of  a  prolate  spheroidal  current  band  in 
free  space  is  shown  when  the  permeability  of  the 
ferromagnetic  prolate  spheroidal  shell  is  allowed 
to  approach  that  of  free  space. 


ADMINISTRATIVE  INFORMATION 

This  work  was  performed  under  Program  Element  11221N,  Project  B0005, 
Task  Area  B0005-SL001,  Work  Unit  2704-110.  The  Project  Director  is 
Mr.  W.  J.  Andahazy,  David  W.  Taylor  Naval  Ship  Research  and  Development 
Center . 


INTRODUCTION 

1* 

As  noted  in  Brown  and  Baker  s  report,  exact  analytical  solutions  of 
Maxwell's  equations  using  classical  formulations  have  been  limited  to  body 
shapes  and  inhomogeneities  that  conform  to  a  few  separable  coordinate 
systems.  Modern  digital  computers  with  large  computational  and  storage 
capabilities  permit  many  electromagnetic  field  problems  to  be  solved  by 
using  a  numerical  solution  to  the  governing  differential  or  integral  equa¬ 
tions  under  a  suitable  choice  of  boundary  conditions.  The  numerical  solu¬ 
tions  of  Maxwell's  equations,  when  used  with  a  complete  description  of  the 
electric  and  magnetic  sources  and  the  constitutive  laws  of  the  media,  can 
be  used  to  describe  completely  the  electric  and  magnetic  fields  produced 
by  the  source,  including  nonsymmetric  geometries,  nonsymmetric  source  dis¬ 
tributions,  and  spatially  varying  media  parameters. 


*A  list  of  references  is  given  on  page  71. 


The  motivation  for  this  work  arose  out  of  the  need  for  solutions  to 
static  ferromagnetic  problems  that  could  be  used  for  comparison  with 
numerical  methods. 


PROLATE  SPHEROIDAL  COORDINATE  SYSTEM 

The  prolate  spheroidal  coordinate  system  can  be  formed  by  rotating 

the  two-dimensional  elliptic  coordinate  system,  whose  traces  in  a  plane 

are  confocal  ellipses  and  hyperbolas,  about  the  major  axis  of  the 
2  3 

ellipses.  ’ 

3 

Flammer  notes  that  it  is  customary  to  make  the  z-axis  the  axis  of 
revolution.  Figure  1  depicts  the  three-dimensional  prolate  spheroidal 
coordinate  system.  In  this  case  the  coordinate  surfaces  are:  prolate 
spheroids  for  n  =  constant;  hyperboloids  of  two  sheets  for  0  =  constant; 
meridian  planes  for  ip  =  constant.  The  prolate  spheroidal  coordinates 
shown  on  Figure  1  are  related  to  rectangular  coordinates  by  the  following 
transformation  equations: 


x  =  a  sinhn  sin0  cosij; 


y  =  a  sinhf|  sin0  sinij/ 


z  =  a  coshn  cos0 


(la) 


(lb) 

(lc) 


where 


o<n<°° 

°<0<Tr 

o<4><2n 

We  have  denoted  the  interfocal  distance  by  2a  and  the  prolate  spheroidal 
coordinates  by  (n,  0,  ip,). 
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PROLATE  SPHEROIDS,  n  -  CONST 


t  =  tt/2 


Figure  1  -  Prolate  Spheroidal  Coordinate  System 
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BASIC  EQUATIONS 

We  can  start  with  Maxwell's  equations  for  classical  magnetostatic 
field  problems 


V  x  H  =  J  (2a)* 

y.  B  =  0  (2b) 

where 

H  magnetic  field  intensity  (A/m) 

_  2 
B  magnetic  flux  density  (T  or  Wb/m  ) 

_  2 
J  electric  current  density  (A/m  )  . 

In  the  general  case  for  ferromagnetic  materials  B  is  a  nonlinear  function 
of  H 


B  =  f (H)  (3) 

wuere  B  is  not  a  single  valued  function  of  H.  The  function  f(H)  depends 
on  the  magnetic  history  of  the  material,  that  is,  how  the  material 
attained  its  magnetization.  This  is  referred  to  as  hysteresis.  It  is 
also  noted  that  any  property  of  a  ferromagnetic  material  has  meaning  only 
if  it  is  considered  together  with  its  complete  magnetic  history. 

In  certain  practical  engineering  problems,  the  variation  in  the  mag¬ 
netic  intensity  is  small,  and  the  functional  relationship  between  B  and  H 
is  approximately  linear.^  For  the  linear  case  where  the  material  is  iso¬ 
tropic,  the  magnetic  induction  B  is  related  to  the  field  intensity  H  by 
the  relationship 


*The  del  operator  V  is  defined  with  respect  to  the  rectangular  coor¬ 
dinate  system  and  is  strictly  valid  in  a  rectangular  coordinate  system 
only.  Very  often  Vx  and  V-  are  used  as  equivalent  symbols  for  curl  and 
divergence  generally.  This  use  is  followed  in  this  report. 
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(4) 


B  =  ^o  (Xm+1)  H  =  MoPrH  =  PH 


where 

magnetic  susceptibility  (dimensionless) 
p  magnetic  permeability  (henry/meter) 

(V1)  =  ux  relative  permeability  (dimensionless) 

free  space  permeability  (4tt  x  10  ^  henry/meter)  . 

This  report  assumes  that  the  ferromagnetic  body  has  isotropic  and  linear 
material  properties.  The  divergenceless  nature  of  the  magnetic  flux  den¬ 
sity  in  conjunction  with  the  fact  that  the  divergence  of  the  curl  of  any 
vector  function  is  zero  allows  the  introduction  of  the  magnetic  vector 
potential  field  (A) 


B  =  V  x  A 


(5) 


where  A  is  the  magnetostatic  vector  potential  function  in  weber/meter.  The 
substitution  of  Equation  (5)  into  Equation  (2a)  gives  the  fundamental 
equation  of  the  vector  potential  of  the  magnetostatic  field. 


-jj  V  x  (VxA)  -  (VxA)  x  V  -jj  * 


(6) 


For  homogeneous  materials,  as  assumed  in  this  report,  the  magnetic 
permeability  is  spatially  invariant.  Hence, 


(7) 


5 


and  Equation  (6)  reduces  to 


V  x  (VxA)  =  yJ 


(8) 


Using  the  vector  identity 


VxVxA=V(V.  A)-  $ A 


(9)* 


Equation  8  becomes 


V(V-  A)  -  0  A  =  yj 


(10) 


The  magnetostatic  vector  potential  is  characterized  by  the  important  prop¬ 
erty  that  its  divergence  can  be  conveniently  chosen  to  be  zero. 


V*  A  =  0 


(11) 


Equation  (10)  reduces  to  the  vector  Poisson's  differential  equation. 


0A  =  -  yJ  (12) 

This  is  the  governing  equation  for  our  calculations. 

The  general  boundary  conditions  to  be  satisfied  at  the  interfaces  of 
stationary  dissimilar  media  may  be  derived  from  the  limiting  integral  forms 
of  Maxwell's  equations  and  are  given  by 


*The  vector  Laplacian  operator  is  designated  by  $  . 
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(13a) 


'12  •  (W  "  0  or  Bnl  =  Bn2 


n12  x  (H2-Hl)  =  Js  or  Ht2  -  Htl  =  Jg  (13b) 


where  the  subscripts  1  and  2  indicate  the  media  under  consideration,  and 
n^2  denotes  the  unit  normal  vector  to  the  interface  and  is  directed  from 
medium  1  into  medium  2.  In  the  case  where  the  materials  are  linear  and 
isotropic  equations  (13a)  and  (13b)  become 


n12  .  (u^-u^)  =  0 


(13c) 


(13d) 


Jg  is  a  true  surface  current  density  that  may  exist  at  the  interface.  At 
an  interface  where  Jg  is  0,  Equations  (13b)  and  (13d)  need  to  be  modified 
accordingly. 


THIN  COIL  SURROUNDING  A  FERROMAGNETIC 
PROLATE  SPHEROIDAL  SHELL 


GENERAL  SOLUTION 

We  now  proceed  to  solve  the  boundary  value  problem  of  a  ferromagnetic 
prolate  spheroidal  shell  of  homogeneous  permeability  y2  surrounded  by  an 
infinitesimally  thin  prolate  spheroidal  current  band  of  constant  current 
density  J.  The  geometry  of  the  problem  suggests  that  a  prolate  spheroidal 
coordinate  system  as  shown  in  Figure  1  be  used  in  the  solution. 

Figure  2,  a  cross  section  of  the  problem  geometry,  identifies  the  four 
regions  of  interest.  The  boundaries  of  the  prolate  spheroidal  shell  are 
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>y  Infinitesimally  Thin  Current  Band 


determined  by  n  =  =  constant  and  n  =  =  constant.  The  dc  current  lies 

in  the  boundary  n  =  constant.  Regions  I,  III,  and  IV  have  a  permeabil¬ 
ity  equal  to  that  of  free  space  UQ»  which  for  convenience  will  be 
labelled  y  .  Ampere's  law  states 


V  x  H  =  J 


(14) 


and  since  B  =  0,  the  induction  B  must  be  the  curl  of  some  vector  field 
A.  The  governing  differential  equation  for  A  when  homogeneous  and  linear 
materials  are  considered  is  from  Equation  (12) 


0A  =  -  UJ  (15) 

The  general  expression  in  prolate  spheroidal  coordinates  for  a  cur¬ 
rent  density  is 


J 


J  e 

n  n 


+  Ve 


+  J  |  e  . 

ip  ip 


(16) 


In  the  problem  presented  herein,  the  current  density  has  only  a  psi  (ip) 
which  means  that  the  vector  potential  has  only  a  psi 


component  jj^(0)e^J 

The  vector  potential  A 


component  A , e , • 

\p  p 


A,e,  is  a  function  of  the  pro- 

ipip- 


lp  r  V  V  I 

late  spheroidal  coordinates  q,  0  li.e.  =  A^(n,8)J  .  The  constant  cur 

rent  density,  which  lies  on  the  boundary  between  regions  III  and  IV,  can 


be  expressed  by  the  function 


0  ,  if  9<0^  or  0>02 

ye>v  if  0i-e-02 


(17) 
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where  Ji/0)  is  e9ual  to  a  constant  J  along  n  =  n3  for  61£0<.02  . 

Therefore,  Equation  (15)  has  only  an  azimuthal  component  and  can  be 
expressed  as 

■  *%<-..>  ■ » *)  «•> 


When  the  psi  component  of  vector  Laplacian  0 is  expanded  in  prolate 
spheroidal  coordinates,  Equation  (18)  can  be  expressed  as  (see  Appendix  A) 


|-  ±  d(sinhn  A^) 

,  _a_ 

l 

a(sin0  A  )■ 

sinhn  an 

a0 

sin0 

d  0 

(in  regions  I-IV) 


(19) 


Applying  the  method  of  separation  of  variables, 
be  expressed  as  the  product  of  two  functions 


let  us  assume  that  A,  can 

'P 


A^  =  H(coshn)  G(cos0)  (20) 

where  H(coshn)  is  a  function  of  coshp  only  and  G(cos0)  is  a  function  of 
cos0  only.  Substituting  this  form  of  the  component  of  the  vector  poten¬ 
tial  A  into  Equation  (19),  we  have  after  separation  of  variables 


+  cothn  ^  -(  p(p+l).  +  — — H  =  0 
dn  dn  \  sinh2n/ 

£|+c°te.d|  +  (p(p+i)  _  1  \G  =  Q 

d0  \  sin  0  / 


(21a) 


(21b) 


10 


where  the  separation  constant  is  p(p+l)  and  p  is  an  integer  from  one  to 
infinity.  It  is  well  known  that  differential  equations  of  the  form 


sV  +  cothn  ai  _  /  p(P+i)  +  -vV' 

dn2  dn  \  sinh2n/ 


=  0 


(22a) 


have  the  general  solution  of  the  form 


H'  =  C-Pm  (coshn)  +  C.,*)”’  (coshn) 
1  p  *  P 


(22b) 


where  and  ate  constants,  and  it  is  known  that  a  differential  equa¬ 
tion  of  the  form 


dV  .a  dG' 
^+C°t6^  + 


P(P+1) 


sin20 


G‘  =  0 


(23a) 


has  the  general  solution  of  the  type 


G'  =  C.Pm  (cos0)  +  C,Qm  (cos0) 
J  p  **  P 


(23b) 


where  and  are  constants.  Pp  and  are  the  associated  Legrendre 

functions  of  the  first  and  second  kind,  respectively.  Comparison  of  Equa- 

2 

tions  (21),  (22),  and  (23)  shows  that  in  Equations  (22)  and  (23),  m  is 
equal  to  1.  This  requires  that  m  always  equal  unity.  The  solutions  of 
Equations  (21a)  and  (21b)  are  expressed  as 


H (coshn)  =  A  P^  (coshn)  +  ®  Qp  (coshn) 


(24a) 
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1  y 


i* 


(24b) 


G(cos0)  =  pj  P1  (cos9)  +  B1'  Q1  (cos0) 
P  P 


The  general  solution  of  Equation  (19)  may  be  formed  from  the  product  of 
solutions  in  Equations  (24a)  and  (24b)  which  yield 


A,  =  H  (coshn)  G  (cos0)  =  E  H  (coshn)  G  (cos0) 

^  p=l  ^ 


(25) 


=  £ 
p=l 


(coshn)  +  B 


(coshn) 


x  A1  P1 
p 


(cos0)  +  B1 


Q  (cos0  )| 


(26) 


For  the  prolate  spheroidal  system,  the  associated  Legendre  functions  of 
the  second  kind  are  infinite  at  cos0  =  ±1  and  as  such  cannot  be  included 
in  a  general  solution  for  a  given  region  which  includes  0  =  0  or  0  =  it. 
Therefore,  in  our  case  the  constant  b'  is  set  equal  to  zero.  Equation  (26) 
reduces  to 


=  E 

p= 


l 


K,  P1 
1  P 


(coshn)  + 


hi 


(coshn) 


P1 

p 


(cos@) 


(27) 


where  ^  and  K2  are  constants  (^  =  AA'  ,  K2  =  AB; )  .  When  the  substitutions 
E,  =  coshn  and  v  *  cos0  are  made  in  Equation  (27)  ,  A^  can  be  expressed  as 


V  E 

P’ 


^  vj®  +  vj®  pj 


(V) 


(28) 


This  is  the  general  form  of  the  psi  component  of  the  vector  potential  that 
will  be  used  to  determine  the  potentials  A^  in  each  region. 
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BOUNDARY  CONDITIONS 

The  form  of  the  component  of  the  vector  potential  in  regions 
(I-IV)  is  determined  from  Equation  (28) .  These  magnetostatic  vector 
potentials  in  regions  I-IV  are: 


A.  T  =  T 

)iil  ^ 

P 


[A  P1(U  P1 

P  P  J  P 


(V) 


>11  =  £ 

p 


>iii 


=  £ 
p= 


;  [vp®  +  cp<®]  pl 

+  Vp1(°]  pp1 


(v) 


(V) 


(29) 


(jiiv 


(V) 


Because  the  potential  must  be  finite  in  each  of  the  regions  I,  II,  and  III 
and  approach  zero  as  (  +  ®  in  region  IV,  the  following  constants  were  set 
equal  to  zero. 

a.  For  A^j  the  constant  associated  with  Q1(5)  P1(v)  was  set  equal  to 


zero  because 


Qp  (O  -*•  00  at  £  =  +1 


b  For  A.  the  constant  associated  with  P1(C)  P.^v)  was  set  equal 


'<J>IV 

to  zero  because  Pp(D  ■+  00  as  (  + 


(We  note  Qp(C)  +  0  as  (  +  ®) 
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The  constants  A  ,  B  ,  C  ,  D  ,  E  ,  and  F  are  to  be  determined  from  the 
P  P  P  P  P  P 

boundary  conditions.  At  each  interface,  the  basic  laws  of  magnetostatics 
(Equations  (2a)  and  (2b))  reduce  to  boundary  conditions  on  B  and  H  that 
can  be  used  to  evaluate  these  six  constants.  The  normal  component  of  B 
across  each  boundary  must  be  continuous,  i.e.,  (B2  -  B^)  '  n^  =  0  where 
the  quantity  n^  is  the  unit  outward  normal  to  the  surface.  This  provides 
the  following  boundary  conditions  which  must  be  satisfied  by  the  solutions 
given  in  Equation  (29)  for  each  region. 


=  BnII  at  n  =  nx 


(30a) 


Bnm  at  n  n3 


(30b) 


Bnm  =  Bmv  at  n  =  n3 


(30c) 


The  eta  or  normal  component  of  the  magnetic  field  (B^)  is  expressed  in 
terms  of  the  vector  potential  as 


_  1  3  <e3A . ) 

Bn  ‘  <V*V,  ‘ 


,(SW) 


(1-v")  A, 


14 


where 


B  =  V  x  A 


2  2 

a(sinh  n  +  sin  0)  (sinhn  sin9) 


2  2  1/2  —  2  2  1/2  - 
e  (sinh  n  +  sin  9)  eQ(sinh  n  +  sin  0)  e ,  sinhri  sin9 

n  t)  ijj 


_d 

an . 


d_ 

dQ 


a_ 

dip 


A ^  sinhn  sin9 


and 


2  2  1/2  1/2 
£  =  coshn,  e^  =  ^  =  a(sinh  n  +  sin  0)  =  a (£  -v  ) 


v  =  cos9,  e^  =  a  sinhn  sin© 


However,  since  the  vector  potentials  in  each  region  are  functions  of  P  (v) 
we  can  simplify  Equation  (30)  to  constraints  on  at  the  interfaces: 


A,  T 

'I'l 

"  \ll 

\ll 

=  Vii 

A 

• 

Vn 

"  \iv 
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(32a) 

(32b) 

(32c) 


'7-rW~ 


.  v  mmrnm '  **s*m*™-~-' - 


The  second  set  of  boundary  conditions  states  that  the  theta  or  tangential 
component  of  H  across  each  boundary  must  satisfy  the  relationship 

ni2  x  =  Jg  (33) 

where  Jg  (which  equals  J^(0)e^)  is  the  real  surface  current  density  in 

the  limit  of  vanishing  width  between  the  two  regions.  Using  the  relationship 

B  =  pH,  Equation  (33)  can  be  expressed  as 


=  ye> 


(34) 


Referring  to  the  curl  in  Equation  (31)  ,  we  can  write  BQ  in  the  form 


Be  =  <VxVe 


(35) 


=  _zL  ?(e3V  =  _ i  r  (r2  _a1/2  A  i 

eie3  .(e-J)1'2  ^  l' 

From  Equations  (34)  and  (35)  the  tangential  components  of  B  in  regions 
I-IV  must  satisfy  the  relationships: 


16 


£ye> 

p 


E 

p 


KGPp>> 


(36c) 


It  is  understood  that  the  summation  runs  from  p  =  1  to  p  =  °°. 

The  general  expressions  for  the  potentials  in  each  region  (Equation  (29)) 
are  then  substituted  into  the  boundary  conditions  (Equations  (32)  and  (36)) 
are  solved  for  the  six  constants  (Ap,  Bp,  Cp,  Dp,  Ep,  and  Fp) .  Since  there 
are  six  equations  with  six  unknowns,  the  potential  in  each  region  can  be 
specified.  The  six  boundary  value  equations  are  presented  below.  The 
index  p  in  the  summation  sign  has  both  even  and  odd  values  and  takes  on 
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values  from  1  to  It  is  noted  at  this  point  that  the  current  density 
J^(0)  must  be  expanded  into  a  set  of  associated  Legendre  functions  in  order 
to  evaluate  the  constants  in  the  vector  potential  (Equation  (29)).  The 
detailed  expansion  is  presented  in  Appendix  B.  The  six  expressions  for  the 
boundary  conditions  are: 


Vp<5i,pJ<u>  -[yX>  +  Vp(Ci)]  pJ<v)  <37a) 


^a(5P-v2)1/2 


(52-i)1/2  (yj«> +  v£<o)pH 


1 


£=?■> 


(37b) 


[bppp(^2>  +  cPQJ(^]  pJ(v)  -  [vJ«. 2> +  y^2>]pJ(v)  (37c) 


(37d) 
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(37e) 


[DpPp<e3>  +  Vj<«3>]  Pp(V)  ‘  V)S>Pp<V) 


p>] 


5-S, 


\1/2,„  a. 


fr)  'a(V2)^  h  [(^  (vJ®  +  ep<co)pp''')] 


(37  f ) 


£-5, 


Jp(0) 


KG  P  (V) 
P  P 


.(cf-v2) 


1/2 


If  we  make  the  following  substitution 


pp«>  ■  3{  (V-*)172  p;«> 


(38) 


Qp(5)  =  ~ 

P  d? 


(c2-i)1/2  qJ<«] 


(39) 


and  perform  simple  algebric  manipulations,  the  six  boundary  conditions 
can  be  simplified  to: 


fWW~ 


v?Ki>  -  yj'v + 


(40a) 
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-uir  ur*j«i 


.  — .f^Li 


(^)[vp<v + «P<‘v]-(y  vx> 


(40b) 


Vp«2>  +  V&V  *  Vp<52>  +  V£«2>  (40c) 


(y  [Vp(52>  +  Vp«2>]  -  (y  [Vp«2>  +  V*S>] 


(40d) 


Vp<53>  +  Ep<(53>  '  Fp<«3> 


(40e) 


(€3) 


(*) 


D  P„(5,)  + 
p  p  3 


EP<(V 


,<9> 


1/2 


pj(v) 


(40f  ) 


where 


Jp(9) 


K  G  p*(v) 
_ P  P 

'(^2)1/2 


The  solution  of  these  six  simultaneous  equations  to  obtain  E  in  terms  of 

P 

known  quantities  gives: 

C  -  ^(MJpn  ppifo)  -  i^K^)  +  p^2> 

^(ww^«2))^(m<«2>)-(~K<«2> 

(41a) 
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where 


(41b) 


M 


0p«2) 

(x1pp<52>  +  qJ«2) 


(41c) 


Jp(0) 


K  yj(v) 


(v=cos0) 


(41d) 


J 


I 

P 


<«3>/ 


(4le) 


jV^,) 

£  P  2 


(4lf ) 


Tlie  numerical  values  for  the  other  five  coefficients  can  be  obtained  from 
the  following  equations: 


VJpI+Epl>]  <42a> 

BP  -  MCP  («b) 
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(42c) 


D 

P 


A 

P 


B  +  C 


P  P 


pJ«i> 


F 

P 


,  +  E 

P 


(42d) 


<42e) 


Since  the  six  coefficients  can  be  determined  for  a  specified  problem 
from  Equations  (41)  and  (42),  the  potentials  anc^  A^iv’ 

in  regions  I  through  IV  can  be  completely  determined.  The  normal  (B^) 
and  tangential  (B^)  to  the  surface  r\  =  constant  (or  £,  =  constant)  com¬ 
ponents  of  the  magnetic  induction  in  each  region  I  through  IV  can  be  deter 
mined  by  using  Equations  (31)  and  (35) . 

THIN  COIL  INTERNAL  TO  A  FERROMAGNETIC 
PROLATE  SPHEROIDAL  SHELL 


GENERAL  SOLUTION 

We  now  proceed  to  solve  the  boundary  value  problem  of  a  ferromagnetic 
spheroidal  shell  of  homogeneous  permeability  u2»  surrounding  an  infinites! 
mally  thin  prolate  spheroidal  current  band  having  a  constant  current  den¬ 
sity  J.  Figure  3  shows  the  cross  section  of  the  problem  geometry.  The 
coordinate  system  shown  previously  in  Figure  1  will  be  used  in  the 
solution.  The  boundaries  of  the  prolate  spheroidal  shell  are  determined 
by  n  =  n3  and  n  =  n2*  The  stea<ty  state  current  lies  in  the  boundary 
n  =  nL  and  between  0jS0£02*  As  in  the  previous  problem,  the  dc  current 
density  has  only  a  psi  component  J  (0)¥^,  and  thus  the  vector  potential 
has  only  a  psi  component  A^e^.  The  vector  potential  is  a  function  of  the 
prolate  spheroidal  coordinates  n  and  0.  The  constant  current  density  is 
expressed  by  Equation  (17)  when  the  boundary  n  is  changed  to  n  =  n^. 
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-  Cross  Section  of  Infinitesimally  Thin  Current 
Surrounded  by  Ferromagnetic  Spheroidal  Shell 


The  governing  partial  differential  equation  has  only  a  psi  component 
and  is  given  by 


=  $A^(n-t8)  =  0  (in  regions  I-IV)  (43) 


When  the  vector  Laplacian  0  A  is  expanded  in  prolate  spheroidal  coor¬ 
dinates,  Equation  (43)  can  be  expressed  as  (see  Appendix  A): 

A  r  1  3  (sinhn  A  .  )  q  i  r  i 

—  - ^  1  +  ^  r^-s-  A-  (sin0  A.)  1  =  0 

an  [sinhn  an  •  J  36  [sine  ae  v  r  J 

Adopting  the  following  notation 


£  =  coshn,  v  =  cos6  (45) 

and  following  the  logic  presented  earlier,  the  solutions  for  A^  in  regions 
I-IV  have  the  general  form 


A  =  L  £kjpJ(5)  +  K2Qp(OjPp(v>  (46) 

BOUNDARY  CONDITIONS 

The  form  of  the  components  of  the  vector  potential  A^  in  each  of  the 
regions  I-IV  is  determined  from  Equation  (46).  These  components  of  the 
vector  potential  in  each  region  are: 


V  ■  £  kpp(#J(v)  <47a> 

p=lL 

Vi  -  Z  [Vi®  +  Vp<0K<v)  <47b> 

p=l 
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>111 


-E[' 

P-il 


yj<« +  v>)|p;<v) 


>]* 


(47c) 


Vv  -£  [v?<5)]p>) 


(47d) 


The  P  functions  are  the  associated  Legendre  functions  of  the  first  kind  of 
P  l 

degree  1  and  order  p,  and  the  functions  are  associated  Legendre  func¬ 
tions  of  the  second  kind. 

At  each  interface  the  basic  laws  of  magnetostatics  reduce  to  boundary 
conditions  on  B  and  H  (see  Equations  (30)  and  (33))  that  can  be  used  to 
determine  related  boundary  conditions  on  A: 


(48a) 


\ll  ~  %III  n  “  n2 


(48b) 


%III  =  %IV  n  “  n3 


(48c) 


£J  (6) 
P  v 


(48d) 
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These  boundary  conditions  are  then  used  to  evaluate  the  constants  in  Equa¬ 
tion  (47) .  Using  Equations  (47)  and  (48)  to  solve  for  the  coefficients 
(where  the  index  p  takes  on  all  values  from  1  to  °°)  we  get: 
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(49c) 


[vJ«2>  +  vJ<v]pJ<'-)  -[ 


Vp«2>  +  V& 


(V) 


("2} 


1  a 


te-2) 


1/2 


&c 


\  1/2 


(e2-1)  (yj«>  +  Vp®)ppm 


c=c 


2  (49d) 


1  \  1  a 


Ul/a(^-v2)1/2  8CL 


(?2-i)1/2  (Vp1(5)  +  \i(») 


5=5, 


[lpPJ<53>  +  Vj'Vttw  -  VJ<E3)PJ<v> 


(49e) 


and  perform  simple  algebric  manipulations,  the  six  boundary  conditions 
reduce  to: 


vX>  -  Vp«i>  +  v>X> 


(52a) 


■(y(  vfa +  v4«i>)  - 


J  (0)a 


f^2): 


pj(v) 


(52b) 


vi(v +  vJ<?2>  ■  spJ«2> + 


(52c) 


(y(Vp«2>  +  Vp<52>)  -(^XVp(52>  +  Vp«2>) 


lpPJ<e3>  +  V>J«3>  ‘  Vp<e3> 


(52e) 


(Si)  Vp<53>  ■  (y(LPPp«3>  +  Vp<S3>)  <52f> 


It  should  be  noted  in  the  above  equations  that  the  current  density 
•J^(9)  was  expanded  into  a  set  of  associated  Legendre  functions  in  order  to 
evaluate  the  constants  in  the  vector  potential  components  (see  Appendix  B) . 
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The  solution  of  these  six  simultaneous  equations  (52a  through  52f)  to 
obtain  in  terms  of  known  quantities  is: 


L  = 
P 


ftt[ 


jPIIpp(v + 


f  57) M  p^(52)  +(yw« 


>p«2> 


(53) 


where 


[u]  . 


~pp(v 

QP«3>  fovU, 


lu2  . 


(54) 


[v]  = 


PJ«2>  +  [“] 


Ft((2> 


(55) 


Jp(6) 


K  G  P  (v) 

_ ELJ?____ 

,uyr 


(56) 


Qp(cl} 

PX) 


(57) 


,11 


-  J 


1 

P  rS> 


(58) 


The  numerical  values  for  the  other  five  coefficients  can  be  obtained  from 
the  following  equations: 


M  =  L  fu"| 
P  PL  J 


Jn+  L  [vl 
P  PL  ^ 


p!;u,) 

N  =  L  -f — —  +  M 
P  P  <£«3>  p 


W*i> 


The  components  of  the  potential  in  regions  I-IV  can  be  determined  since 

the  coefficients  H  ,  I  ,  K  ,  L  ,  M  ,  and  N  can  be  calculated  for  a  speci- 
P  P  P  P  P  P 

fic  problem.  The  normal  (B  )  and  tangential  (BQ)  components  (to  the  sur- 

n  o 

face  n  =  constant  or  (  *  constant)  of  the  magnetic  induction  in  each  region 
1-IV  can  be  determined  by  using  Equations  (31)  and  (35)  . 
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APPENDIX  A 


DERIVATION  OF  THE  VECTOR  LAP  LAC  IAN  0*^,0) 

We  note  that  a  distinction  is  made  between  the  Laplacian  operating  on 

2  ^ _ 

a  scalar  V  tj>  and  the  Laplacian  operating  on  a  vector  ^A(n,0)  [i.e., 

2  _ 

V  A(n,9)] .  The  vector  Poisson's  equation  in  rectangular  coordinates  can 

2 

be  treated  as  three  uncoupled  scalar  equations  where  V  A^  =  for 
i  =  x,y,z.  However,  if  the  vector  Poisson's  equation  is  resolved  into 
orthogonal  components  in  other  coordinate  systems  the  differential  opera¬ 
tion  mixes  the  components  together  giving  coupled  equations. 

The  vector  Laplacian  ^A  can  be  derived  by  using  the  well-known 
vector  identity 

V  x  (VxA)  =  V(V-A)  -  $A  (A-l) 

In  the  case  herein  the  Coulomb  gage  was  chosen  (V*A=0);  therefore,  the  fol¬ 
lowing  identity  applied  in  our  magnetostatic  work 


^A  =  -  V  x  (VxA) 


(A-2) 


The  expression  in  prolate  spheroidal  coordinates  for  0A^(n,0)  will  now  be 
derived  since  only  the  psi  component  of  the  vector  potential  A  (ri*0)  exists 
in  this  problem. 

Taking  the  curl  of  A  results  in  the  expression 

v  x  A^(n,0)  = 


e , sinhn  sin0 

A  sinhn  sin0 
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- ^  -p-  sinhn  A, 

'tan  sinhn  dn  4 


(A-7) 


(sinh  n+sin  Q )  —  /  d 

a2 (sinhn  sin0) 


+  ^0  (sin6 


Now,  if  $A^  is  set  equal  to  zero,  then  the  vector  potential  A^(n»0)  will 
satisfy  the  equation  (see  Equation  (19)  in  text  of  report) . 


^  ^(sinhn  V ]+  ^  [sine  ^(sin0  v]“  0 


(A-8) 


It  is  also  interesting  to  note  at  this  point  that  V  ’A^n,®)  =  0  as  it 
must,  since  the  Coulomb  gage  was  chosen.  The  divergence  of  the  vector 
A  in  prolate  spheroidal  coordinates  is 


V  •  A 


a (sinh 


{1/2 

— - —  |7sinh2rr+sin20)  sinhn  A  1 

sinhn  dn  L'  J 


(A-9) 


¥I¥0  ^  [(si«hVsin20)  sin©  A0, 

t 


sinn  sin0  gij; 


Since  the  two  components  A  and  AQ  are  zero  in  our  case,  and  A  is  a 

n  _  v 

function  of  n  and  0  only,  then  V-A(jJ(n,0)  is  equal  to  zero  as  required. 
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APPENDIX  B 


When  J.(0)  is  set  conveniently  to  a  constant  J,  G  can  be  expressed  as 

ip  P 


=  (2p+l)a 
P  2p(p+l) 


yrm  n  1/2 

^sinh^n^+sin^o)  Pp(cos0)sin0d0 

r\ 


(B-3) 


Since  the  current  in  this  problem  extends  from  0^  to  0^  the  expression  for 

G  can  be  further  modified  to  read 
P 


f  ^sinh2ri1+sin^0^  (cos0)  sin0d0 


(B-4) 


By  making  the  following  substitutions 


£  =  coshn  ,  V  =  cos0 


(B-5) 


G  can  be  conveniently  written  as 
P 


G  =  -  (2P+l)a 

P  2p(p+l) 


rv  2  ! 

f  (& 


P  (v)dv 
P 


(B-6) 


This  integral  can  be  evaluated  with  much  difficulty  in  closed  mathematical 

4 

form.  Purczynski  evaluated  these  integrals  for  C£p  (even  terms)  and 
c2p+1  (°dd  terms)  in  terms  of  very  complicated  series  expansions. 

In  the  computer  work  to  be  performed,  the  integrals  for  Gp  will  be 
evaluated  numerically,  the  reason  being  it  has  been  found  that  a  large 
number  of  coefficients  were  necessary  to  obtain  a  good  fit  for  the  current 
expanded  in  terms  of  associated  Legendre  functions  for  the  analogous 
spherical  shell  problem. 
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APPENDIX  C 


DERIVATION  OF  THE  COEFFICIENTS  FOR  THE  VECTOR  . 

POTENTIALS  FOR  A  THIN  COIL  SURROUNDING 
A  FERROMAGNETIC  SPHEROIDAL  SHELL 

In  this  appendix  the  coefficients  are  derived  for  the  vector  potential 
in  regions  I-IV  for  a  ferromagnetic  spheroidal  shell  surrounded  by  an 
infinitely  thin  current  band.  For  a  detailed  discussion  of  the  ferromag¬ 
netic  problem  see  the  section  in  the  text  of  this  report  entitled  "Thin 
Coil  Surrounding  a  Ferromagnetic  Prolate  Spheroidal  Shell."  The  magnetic 
vector  potentials  in  each  region  are  given  by: 


V  ■'j'1[Vp<£,]pp(v) 

GO  ^ 

Vi  •  £  [yj«>  +  7)P(5)]PP<V) 

P=1  u 

Vi  ■  £  [■ yj«>  -  Epc*p(«>]pp<v) 

P=1 


(C-la) 


(C-lb) 


(C-lc) 


(C-ld) 


The  coefficients  in  Equations  (C-la)  to  (C-ld)  are  obtained  by  substituting 
these  equations  into  Equations  (C-2a)  to  C-2f ) . 


\l  =  A\pII 

A\pII  ~  \ III 

A  s  A 

>111  >IV 


at  r) 


at  n 


at  n 


(C-2a) 


(C-2b) 


(C-2c) 
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-(sr)v^3>  *  (yvte> + (y  v$«3> 

1  P  (v) 

p 


1/2 


(C-3f ) 


where 


k  G  P  (v) 

Jp '  .(c2-J)1/2 


These  algebraic  equations  provide  six  simultaneous  equations  with  six 

unknowns,  and  they  can  be  solved  for  the  coefficients  A  ,  B  ,  C  ,  D  ,  E  , 

P  P  P  P  P 

and  F  . 

P 

After  solving  Equation  (C-3a)  for  A^,  we  obtain 


A 

P 


B 

P 


+  C 

P 


(C-4) 


The  solution  of  A  from  Equation 


(C-3b) 


is 


A  = 
P 


C 

P 


l,2\p5«i)J 


(C-5) 


By  equating  Equations  (C-4)  and  (C-5)  and  solving  for  B 


P’ 


one  obtains 


C 


P 


.W2\Pp^!)/ 

P1^) 
_ E _ L- 

B  = 
P 


where 


Solving  Equation  (C-3e)  for  F 


we  derive  the  result 


F 

P 


D  PP1<;3)  ,  c 
"  QJ(53)  ” 


(C-8) 


(C-10) 
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where 


/  2  2  \'L^ 
-Vp  a(q-V  ) 

PP(^3> 

pJw»qJ«3) 

/ 

qJ«3> 

<£«3> 

- 

Solving  Equation  (C-3c)  for  Cp  and  using  Equations  (C-7)  and 
(C-10),  we  arrive  at  the  result 


EPy<.2> 


(WPp«2>  +  <£«2>)  (MpJ«2>  +  «J«2>) 


This  can  be  expressed  as 


where 


C  =  J11  +  E  [  z] 
P  P  p  L  J 


II 


Jp?P(t2> 


[WPJ«2>  + 


(C-ll) 


(C-12) 


[WpJ«2)  +  oJ«2)] 
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'^'T'  tP - 


Using  Equation  (C-3d)  to  solve  for  and  substituting  Equations  (C-7), 
(C-10) ,  and  (C-12) f  we  arrive  at  the  expression 


E 

P 


Pp^2> 


-W'“) 


«p«2> 


jV(l) 

p  p  2 


Pp«2> 


)  +  k  (W  ^ 


W  <£M  '  TT  >ip<«2> 


(C-13) 


The  constants  are  thus  determined  by  using  Equations  (C-4) ,  (C-7),  (C-8) , 
(C-10),  (C-12),  and  (C-13).  For  convenience,  (C-4)  ,  (C-7),  (C-8), 

(C-10) ,  and  (C-12)  are  listed  below 


C 

P 


jI1  +  E  H 

p  pL  J 


B 

P 


D 

P 


J 


I 

P 


A 

P 


B 

P 


+  C 

P 


pX> 


(C-14) 


F 

P 


D 

P 


2> 


+  E 

P 
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APPENDIX  D 


DERIVATION  OF  THE  COEFFICIENTS  FOR  THE  VECTOR 
POTENTIAL  FOR  A  THIN  COIL  INTERNAL  TO 
A  FERROMAGNETIC  SPHEROIDAL  SHELL 

In  this  appendix  the  coefficients  are  derived  for  the  vector  poten¬ 
tials  in  regions  I-IV  for  an  infinitely  thin  current  band  surrounded  by  a 
ferromagnetic  spheroidal  shell.  For  a  detailed  discussion  of  the  ferro¬ 
magnetic  problem  see  the  section  in  the  text  of  this  report  entitled  "Thin 
Coil  Internal  to  a  Ferromagnetic  Prolate  Spheroidal  Shell."  The  magnetic 
vector  potential  in  each  region  is  given  by: 


V  ■  ?1[hppp<55pp<v)  <D'la> 

V  ■  |[vJ<?) +  kpqp(5,]pJ(v)  (D‘lb> 

Vi  ’  fjvj®  +  mpqJ«>]pp<u)  <D-lc> 

00  r 

Vv  ■  fi[Vp®]vv)  <D"ld) 


The  coefficients  in  Equations  (D-la)  to  (D-ld)  are  obtained  by  substituting 
these  equations  into  the  boundary  conditions  (Equations  (D-2a)  to  (D-2f) . 


V  =  Vi  at  *  = 
Vi  =  Vii  at  n  “  n2 

Vn  "  Vv  at  n  "  n3 


(D-2a) 


(D-2b) 


(D-2c) 


45 


(D-2f) 
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After  appropriate  substitution  of  Equations  (D-la)  to  (D-ld)  into  Equations 
(D-2a)  to  (D-2f ) ,  the  following  boundary  value  equations  are  obtained. 


vi«i>  ■  yJ<{i> +  vJ<5x) 


(D-3a) 


1/2 

-fc)  (Vp«  > + 


(D-3b) 


IpPp<«2>  +  KpQP<52>  '  LpPJ«2>  +  »pQp«2> 


(D-3c) 


(u2)  (“pW  +  MpQp<52>)  ■  (y  (IpPp<S2)  +  ltpl)p<?2))  (D-3d) 


Vj'V  +  Vj«3)  -  NpqJ({3) 


(U-3e) 


fc)  Vp'V  ‘  (^PPP«3>  +  Vp«3>) 


(U-3f ) 


These  algebraic  equations  provide  six  simultaneous  equations  with  six 

unknowns,  and  they  can  be  solved  for  the  coefficients  H  ,  I  ,  K  ,  M  ,  and 
N  .  P  P  P’  p’ 

P 

Solving  Equation  (D-3a)  for  Hp  in  terras  of  I?  and  Kp,  we  arrive  at 
the  equation 


H 

P 


I  +  K 


P  P 


<«1> 

p;<v 


(D-4) 
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Solving  Equation  (D-3b) 


for  H  , 
P 


we  obtain  the  equation 


H 

P 


I  +  K 
P  P 


WpP1(v) 


(D-5) 


Equating  Equations  (D-4)  and  (D-5)  and  solving  for  K  yields 


K 

P 


(D-6) 


where 


/  2  2\^2-| 
Vpa(q-v ) 

/ 

’«i«i> 

i^i> 

/ 

Now  solving  Equation  (D-3e)  for  N  ,  we  obtain 

P 


N 

P 


L 

P 


PS«3> 

^«3> 


+  M 

P 


(D-2) 


Additionally,  by  solving  Equation  (D-3f)  for  N 


one  obtains  the  expression 


N 

P 


yX> +  Vp«3>) 


«p«3> 


(D-8) 
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liy  equating  Equations  (D-7)  and  (D-8)  ,  we  may  obtain  in  terms  of  Lp 


M  =  L  [ul 

P  pl  J 


(D-9) 


where 


Pp«3>  Pp«3>  lVl\ 


Sj«3>  <<«3 


)  ' 


K  J 


Solving  (D-3c)  for  Ip  and  using  Equations  (D-6)  and  (D-9),  we  derive  the 
result 


I  =  +  L  [v] 

P  P  Pl  J 


(D-10) 


wliere 


jii , .  ji  fei 

P  P  „!/■£-  % 


P  PJ«2> 


pl«2>  +  [“]<!P(52) 


PJ«2> 


Now  if  we  use  Equation  (D-3d)  and  substitute  from  Equations  (D-6),  (D-9), 

and  (D-10)  the  constant  L  can  be  determined  to  be 

P 


(^)(Jplpp«2>  +  JX«2>1 
(il)Mpp<52>  +  (^)pp«2)  +(^)W« 


(D-ll) 
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The  constants  have  now  been  found.  After  the  numerical  value  for  L  is 

P 

calculated  on  the  computer  for  a  specific  problem,  the  numerical  values 
for  the  other  coefficients  can  be  obtained  from  the  following  equations: 


K 

P 


J 


I 


P 


N 

P 


PX> 

L  -? — —  +  M 

P  0p«3>  P 


H 

P 


I  +  K 


P  P 


pJ«i> 


(D-12) 
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APPENDIX  E 


DETERMINATION  OF  THE  MAGNETIC  VECTOR 
POTENTIAL  FOR  AN  INFINITESIMALLY  THIN 
PROLATE  SPHEROIDAL  CURRENT  BAND 

In  this  appendix  the  Potentials  A^  in  the  inner  region  and  A^^.  in 
outer  region  are  derived  for  the  infinitesimally  thin  current  band  in  a 
homogeneous  medium  of  permeability  (see  Figure  E-l) . 

The  potential  in  the  inner  region  A^j.  and  the  outer  region  A^j.  of 
the  infinitesimally  thin  current  band  problem  are  solutions  to  the  vector 
Laplace's  equation  $A  =  0.  These  solutions  can  be  expressed  as: 


A,t  =  E  A  P1(5)  P1  v 

ipi  ^  p  p^7  pv 
p=l 


(E-la) 


V  ‘  ^  plP  M 


(E-lb) 


The  coefficients  A^  and  F^  are  determined  from  the  boundary  conditions  of 
the  problem.  After  algebraic  manipulation  such  as  with  Equations  (32)  and 
(36)  in  the  text,  the  boundary  conditions  for  the  normal  component  of  B 
and  the  tangential  component  of  H  becomes: 


A\pl“  A\JjII  at  n  =  n! 


(E-2a) 


-(H  1  .1/2  —  r(e2-^)1/2A»xxl 

'".(ej-v2)  »«L  1 


(k)  a(5/v2)172  k 


(E-2b) 


Jp(6) 
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Using  the  relationship  £  =  coshh  and  v  =  cos0  and  substituting  the  expres¬ 
sions  for  A^j.  and  A^j.  (Equations  (E-la)  and  (E-lb))  into  the  boundary 
value  equations  (Equations  (E-2a)  and  (E-2b))  provides  us  with  the  follow¬ 
ing  algebraic  equations  for  the  coefficients: 


(E-3a) 


The  following  substitutions  were  used  to  simplify  the  above  expressions. 

pp«>  •  ^  [(52-l)  1/2pJ«>J  <p-‘) 

Qp(5)  -  ^(52-i)1/2qJ(S)1  (e-5) 


These  equations  are  solved  for  Ap  and  Fp  by  simple  algebraic  manipulation. 


A 

P 


f 

P  P^l> 


(E-6a) 


1  2 

piVe>  *(«!-''  ) 

l\i «i>  <«i>] 

.  pp«i>pJ<'’>  -1/ 

LpJ<«i>  pp«i>J 

(E-6b) 
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I 


The  potential  and  A^j.  are  determined  by  substituting  the  expression 
for  Ap  and  into  Equations  (E-la)  and  (E-lb) 


APPENDIX  F 


REDUCTION  OF  THE  MAGNETIC  VECTOR  POTENTIAL 
FOR  A  THIN  COIL  SURROUNDING  A  FERROMAGNETIC 
SPHEROIDAL  SHELL  TO  THAT  OF  A  THIN  COIL  IN 
FREE  SPACE  WHEN  IN  THE  LIMIT  y2  EQUALS  U  x 

In  this  appendix,  the  coefficients  A  ,  B  ,  C  ,  D  ,  E  ,  and  F  for  the 

P  P  P  P  P  P 

potentials  are  evaluated  for  the  system  consisting  of  a  ferromagnetic  shell 
with  permeability  y2  surrounded  by  an  infinitesimally  thin  current  band 
in  a  homogeneous  medium  with  permeability  y^  in  the  limit  as  y2  = 

These  coefficients  are  utilized  in  Equation  (29)  in  the  section  of  this 
report  entitled  "Thin  Coil  Surrounding  a  Ferromagnetic  Prolate  Spheroidal 
Shell."  The  variables  are  defined  in  Figure  2  located  in  the  text.  When 
y2  us  set  equal  to  y  the  problem  reduces  to  that  of  finding  the  potentials 
in  the  two  regions  of  a  simple  current  band  (see  Figure  1-E  in  Appendix 
E  ,  since  the  ferromagnetic  shell  will  now  have  a  permeability  y^  equal 
to  that  of  the  homogeneous  medium  with  permeability  y^. 

In  this  limit  the  coefficients  should  assume  the  following  form: 


A  =  B  =  D 
P  P  P 


(F-la) 


C  =  E  =  0 
P  P 


(F-lb) 


and  where  A  and  F  should  reduce  to  the  coefficients  for  the  potentials 
P  P 

in  the  two  regions  for  the  spheroidal  band  problem  (see  Appendix  E) .  If 
the  coefficients  assume  this  mathematical  form  it  will  prove  that  the 
mathematical  form  of  the  coefficients  for  the  spheroidal  shell  surrounded 
by  a  thin  current  band  are  mathematically  correct. 

The  mathematical  solution  for  E^  in  terms  of  known  quantities  was 
derived  in  Appendix  C  and  was  reported  in  the  text  of  this  report 
(see  Equation  (41a)). 


(F-2a) 


E 

P 


-  -  k  ^S>  ♦  ^ 


(F-2b) 


M  = 


i«2> 


[*]PJ«2>  +  «X> 


pvs2' 


(F-2c) 


J1- 

P 


U, J  (0)a 
_  1  P 


M 


1/2 


pJ(v)qA(^) 


rif£.!^V 

qJ(53)  q£(S3> 


(F-2d) 


.II 


Jpfp«2> 


WpJ(52)  +  qjcc,) 


p  ^2 


(F-2e) 


The  coefficient  Ep  will  now  be  evaluated  when  the  limit  is  taken  with 
M2  =  yx  which  cause  [x]  to  approach  infinity  («) .  Also  the  expression 

for  JpX  is  substituted  into  Equation  (F-2a) . 


56 


limit 


The  expression  for  B 

P 


(see  Equations  (42b)  and  (41b))  is 


BP  ■  w  cP  -  +  epmm 


where 


C 

P 


The  expression  for  B  when 

P 


equals 


can  be  expressed  as 


- 

B 

P 

..  . 

=  limit 

Vyi 

M 

E 

P 

U2 

=U1 

. 

=V 

[zj  j+  limit 


B 

P 


J 


I 

P 


1 


where 


limit 

V^i 


[mw]‘  <«2)  z> 


and 


limit 


M2-nj 


1 


rrxijnl=  iimit 

1 L  J  p 

J  [x)->  °° 

(F-4a) 


(F-4b) 


(F— 5 ) 
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and 
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The  expression  for  A 

P 


is  (see  Equation  (42d)  in  the  text) 


A 

P 


+  C 


^«i> 

P 


(F-8) 


The  expressions  for  A^ 


when  y2m 


is 


A  I  =  B 

+ 

C 

P  P 

P 

M2=yi  y2 

=yi 

.y2 

=yl.. 

(F-9) 


where 


P 

Vyi 


=  0 


y2=yi 


V2=Ul 


The  expression  for  F 

P 


(see  Equation 


(42e)  in  the  text)  is 


F 

P 


V>g3> 

«J«3> 


+  E 


(F-10) 


The  expression  for  F  when 

P 


M2=Ml 


is 


P2-Ml 


y2-yl 


PJ^3> 

— —  +  E 

P 


Vyi 


60 


Vwi 


2  2  \1/2 
V,  a(Vv  ) 

/ 

«J«3> 

QP«3> 

Fpl53>pJ<v) 

/ 

*<3>J 

(F-ll) 


where  E 


=  0  from  Equation  (F-3) 


U2=Pi 


This  means  that  in  the  four  regions,  the  potential  used  in  Equation  (29) 
of  the  main  text  of  this  report  reduce  when  to  ^orm 


>1’  II’  III 


=  E 


p=i 


pJ(5)pJ<v) 


P2-Pl 


(F-12a) 


%IV  ^ 
p=l 


F 


P 

VM1 


y5>Pp<v> 


These  are  the  solutions  for  the  potentials  of  the  current  band  in  a  region 

of  space  with  homogeneous  permeability  (see  Equations  (E-la)  and 

(E-2a)).  We  now  have  the  solution  for  the  two  potentials  in  regions  I 

and  II  (A^p  and  A^.^,  respectively)  for  the  simple  current  band 

problem.  This  indicates  that  the  form  of  the  coefficients  A  through  F 

P  P 

is  correct. 

The  mathematical  expressions  for  A  and  F  which  has  been  evaluated 

P  P 

in  the  limit  as  y2=y 2»  will  be  compared  with  the  coefficients  Ap  and  Fp, 
respectively,  for  the  two  regions  of  the  current  band  problem  (see  Appen¬ 
dix  E)  . 

Ap  from  Equation  (F-9)  is 


(F-13) 
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A  from  Equation  (E-6a)  which  we  will  denote  here  by  A'  is 

P 


A' 

P 


piJP  ‘kW', 

1/21 

1/ 

[  pp«l  >p>> 

PJ«1>. 

,pJ«i>  #{i>. 

(F-14) 


By  simple  algebraic  manipulation,  the  identity 


l 

QP(5i)l 

l 

PP<«1>  pJ«i> 

<£«i> 

.pJ«i> 

PX> 

pS> 

PJ«1>. 

.qp<£i> 

(F-15) 


can  be  changed  to  read 


r%1(v  foy 

LPJ«1>  Pp«l> 


'fel.  fe' 

."p«i>J  Lpp«ii 


rpp«i>  pX>~1 


(F-16) 


substituting  Equation  (F-16)  into  (F-14)  we  have 


A* 

P 


-Vp4i-"2j 

<»p<el>pp«v» 


\l/2 


rp<H>  PP«X> 

L<«1>  <«i> 


(F-17) 


comparing  Equation  (F-17)  to  (41e)  shows  that 
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( F— 18 ) 


Now  F  from  Equation  (E-6b)  which  will  be  denoted  by  F7  is 
P  P 


F7 

P 


,  ,  1/2  . 
v,  ) 

7 

«J«1>  <<«!> ' 

/ 

,Pp«l>  Pp«!)  . 

which  is  identical  to  F  as  defined  in  Equation  (F-ll)  so  that  F  =  F 7 

P  P  P 


Thus,  the  mathematical  expressions  for  A  and  F  (Equations  (F-13)  and 

P  P 

(F-ll) ,  respectively)  for  the  ferromagnetic  spheroidal  shell  surrounded 
by  a  thin  current  band  in  the  limit  as  nte  the  same  as  the  coeffi¬ 

cients  Ap  and  Fp  (see  Equations  (E-6a)  and  (E-6b),  respectively)  for 
the  vector  potentials  in  the  regions  of  the  current  band  in  free  space 
(see  Appendix  E) . 

It  is  noted  that  when  making  the  comparison,  £  must  be  set  equal  to 
£].•  I'or  comParis°n,  the  coefficients  for  the  current  band  problem  are 


A 

P 


F 

P 


p£«i>  p 


(F-19) 


F 

P 


r  2  i\  ~ 

VP  “(h-"  ) 

/ 

pp«1)pJ<',> 

/ 

,PJ«1>  Pp«l>. 

(F-20) 


and  the  coefficients  for  the  ferromagnetic  shell  problem  with  p  =p  are 


F  VS3} 
P  - 

Pp«3J 


(F-21) 
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APPENDIX  G 


REDUCTION  OF  THE  MAGNETIC  VECTOR  POTENTIAL  FOR  A  THIN 
COIL  INTERNAL  TO  A  FERROMAGNETIC  PROLATE  SPHEROIDAL 
SHELL  TO  THAT  OF  A  THIN  COIL  IN  FREE  SPACE  WHEN 
IN  THE  LIMIT  y  EQUALS  y  1 

In  this  appendix,  the  coefficients  H  ,  I  ,  K  ,  L  ,  M  ,  and  N  for 

P  P  P  P  P  P 

the  potentials  are  evaluated  for  the  system  consisting  of  an  infinite¬ 
simally  thin  current  band  surrounded  by  a  ferromagnetic  shell  with 
permeability  y ^  in  a  homogeneous  medium  with  permeability  y^  in  the  limit 
as  y2=y^.  These  coefficients  are  utilized  in  Equation  (47)  in  the  section 
of  this  report  entitled  "Thin  Coil  Internal  to  a  Ferromagnetic  Prolate 
Spheroidal  Shell."  The  variables  are  defined  in  Figure  3  located  on  page 
23  of  this  report.  When  is  set  equal  to  y^  the  problem  reduces  to  that 
of  finding  the  potential  in  the  two  regions  of  a  simple  current  band 
(see  Figure  1-E  in  Appendix  E) ,  since  the  ferromagnetic  shell  will  now 
have  a  permeability  y^  equal  to  that  of  the  homogeneous  medium  with 
permeability  y^. 

In  this  limit  the  coefficients  should  assume  the  following  form 


K  =  M  =  N  (G-la) 

P  P  P 


L  =  0 
P 


(G-lb) 


and  where  H  and  N  should  reduce  to  the  coefficients  for  the  potentials 
P  P 

in  the  two  regions  for  the  spheroidal  band  problem  (see  Appendix  E) .  If 
the  coefficients  assume  this  mathematical  form  it  will  prove  that  the 
coefficients  for  the  current  band  surrounded  by  a  spheroidal  shell  are 
correct. 

The  mathematical  solution  for  L^  in  terms  of  known  quantities  was 
derived  in  Appendix  D  and  was  reported  in  the  text  of  this  report  (see 
Equation  (53)). 


(y[jpipp«2> + jx«2> 


(G-2a) 
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where 


[u]  = 


Pp«3> 

<«3>/M 

^«3> 

r^-ii 

l'w2  J 


(G-2b) 


PJ«2> 


(G-2c) 


k  Gr,K^ 
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(G-2d) 


JX  = 
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/- 2  2\1/21 
1,iJ»  nerv ) 

/ 
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/ 

Pp<^l). 

(G-2e) 


J11  =  -  J 


1  QP(^2) 

P  Pp(V 


(G-2f ) 


The  coefficient  L  will  now  be  evaluated  when  the  limit  is  taken  with 
P 


^2=^1  which  causes  [u]  and  [v]  to  approach  infinity  (°°) 


V^i 


=  limit  < 
u 

Vl-*-  00 


1 

^1 

.J”PP«2 

+ 

JdQd(52) 

p  p  2 

’p«2>  + 
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Pp(52J 

(G-3) 
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(G-12) 
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y2“Ml 
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I 
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It  is  noted  that  this  expression  is  identical  to  F  of  Equation  (E-6b) . 
The  coefficient  Np  is  expressed  as  (see  Equation  62) 


N 

P 


-L  ^  ,  „ 
p  <«3>  " 


(G-13) 


and  since  L  was  shown  to  be  zero  in  Equation  (G-3) ,  N  is  equal  to  M  . 

‘  P  P 

The  expression  for  Kp  is  defined  in  Equation  (59)  to  be 


K 


P 


(G-14) 


which  is  equal  to  N  and  M  . 

P  P 

The  last  coefficient  is  Hp  which  can  be  expressed  as  (see  Equation  63) 


H 
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I 
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+  K 

P 


PJ«1> 


(G-14) 


12  2 
Vp(9)  a(q-v ) 

^(«i>' 

/ 

^1 

/ 

pJ«i>  pp«i> 

- 

(G-16) 


This  is  identical  with  the  expression  for  A  in  Equation  (E-6a)  of 
Appendix  E. 
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Thus,  the  mathematical  expressions  for  and  (Equations  (G-16)  and 
(G-13) ,  respectively)  for  the  thin  current  band  surrounded  by  a  spheroidal 
shell  in  the  limit  as  y  =u,  are  the  same  as  the  coefficients  A  and  F 
(see  Equations  (E-6a)  and  (E-6b) ,  respectively)  for  the  vector  potential 
in  the  regions  of  the  current  band  in  free  space  (see  Appendix  E) . 

For  comparison,  the  coefficients  for  the  ferromagnetic  shell  problem 
with  y2=^i  are 


and  the  coefficients  for  the  current  band  problem  are 
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